One-dimensional Hubbard model at quarter filling on periodic potentials 
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Using the Hubbard chain at quarter filling as a model system, we study the ground state properties 
of highly doped antiferromagnets. In particular, the Hubbard chain at quarter filling is unstable 
against and 4fcF-periodic potentials, leading to a large variety of charge and spin ordered ground 
states. Employing the density matrix renormalization group method, we compare the energy gain of 
the ground state induced by different periodic potentials. For interacting systems the lowest energy 
is found for a 2fcF-periodic magnetic field, resulting in a band insulator with spin gap. For strong 
interaction, the 4fe_F-periodic potential leads to a half- filled Heisenberg chain and thus to a Mott 
insulating state without spin gap. This ground state is more stable than the band insulating state 
caused by any non-magnetic 2fc_F-periodic potential. Adding more electrons, a cluster-like ordering 
is preferred. 
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HOLE DOPED ANTIFERROMAGNETS 



The study of low-dimensional doped antiferromagnets 
was triggered by investigations of the phase diagram of 
high temperature superconductors. Upon doping, they 
show a transition from an antiferromagnetic to a super- 
conducting state. On the other hand, the investigation 
of one-dimensional compounds is also an active held. A 
simple form of low doping in a Cu-chain is the replace- 
ment of a few Cu ions by non-magnetic ions like Zn. In 
this case, the system is well described by the Heisen- 
berg model with some spins removed The charge 
and spin order of highly doped compounds, on the con- 
trary, are not understood by far. A prototypical ma- 
terial revealing strongly doping dependent charge order 
and antiferromagnetism in a less than half-filled band is 
the spin chain compound (Ca,Sr)i4Cu2404i. This com- 
posite crystal contains two different structural compo- 
nents along the c-axis: CU2O3 ladders and CuC>2 chains. 
Whereas the ladders contain only Cu 2+ -ions, nearly all 
Cu 3+ -ions are located on the chains 0. Substitution of 
calcium on the strontium sites leads to a transfer of holes 
from the chains into the ladders. For the Ca-rich com- 
pounds one generally assumes one hole per formula unit 
on the ladders and five on the chains. Near the quarter- 
filled band antiferromagnetic order is established [|[. 

DFT-LDA calculations 0], using the ASW method @], 
reveal that the ladder and chain substructures can be 
treated individually. In particular, the density of states 
of the chains shows a partially filled valence band, as 
illustrated in Fig [TJ The origin of antiferromagnetism 
in these compounds is still unclear, and the doping de- 
pendence of the charge order and spin gap is discussed 
controversially @. Recent studies on the basis of the 
Heisenberg model show that a simple spin model is not 
appropriate to describe the magnetization [3] , indicating 
the necessity of studying more advanced models, as e. g. 
the Hubbard model. 

In this article we address the stability of the pro- 
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FIG. 1: (Color online) Partial Cu 3d density of states for the 
Cu02-chains of Sri4Cu2404i. 



posed charge order (cluster-like versus staggered) Q in 
(Ca,Sr)i4Cu2404i, considering the competition between 
kinetic energy, interaction and periodic distortions, on 
the basis of the Hubbard model. Since no periodic lat- 
tice distortion - neither with two-site nor with four-site 
periodicity - is found in the crystal structure [§] we con- 
sider only periodic potentials. Nevertheless, our results 
are relevant for all correlated materials with quarter filled 
bands. Thus we aim at understanding a whole class of 
materials on a common basis. 
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II. PERIODIC DISTORTIONS IN 
INTERACTING CHAINS 

Charge order - independent of magnetic ordering - 
in less than half-filled bands is found, for example, in 
organic charge transfer salts, like the Bechgaard salts 
(TMTSF) 2 X, with X=PF 6 , AsF 6 , C10 4 , Re0 4 , or Br, or 
their sulfur analogs (TMTTF) 2 X. These systems, which 
consist of stacks of organic molecules forming weakly cou- 
pled one-dimensional chains, exhibit a large variety of low 
temperature phases [l(j, EH ■ They have been studied us- 
ing a purely electronic model, namely the extended Hub- 
bard model with next-nearest-neighbor interaction. Re- 
cent results have been obtained by mean-field approaches 
[l2l ] and exact diagonalization [l3j]. On the other hand, 
charge order is often connected to a structural transition 
and hence to clcctron-phonon coupling. The coupling of 
a one-dimensional metal to an elastic lattice results in an 
instability towards a periodic lattice distortion commen- 
surate to the band-filling, known as Peierls transition. 
Within a lattice model, the Peierls transition is captured 
by a modulated hopping term in the Hamiltonian of the 
non-interacting system [Tj] . A periodic potential, instead 
of a modulated hopping, has similar effects on the elec- 
tronic structure of the one-dimensional metal, leading to 
a band insulator. 

The interplay of electron-lattice and electron-electron 
interaction is conveniently studied in the framework of 
the Hubbard model. With interaction, the Peierls and 
the so-called ionic model show different ground-state 
phase diagrams [HI]. Many results are available for the 
half-filled Hubbard model with periodic hopping or po- 
tential, for which the transition from a band insulator to 
a Mott insulator has been studied in [U HH H3, [H, H . 
On the other hand, away from half filling no insulator- 
insulator-transition is present, but a variety of ordering 
processes have to be taken into account [2a]. The sta- 
bility of the different ground states of the quarter-filled 
Hubbard model with respect to a periodic potential has 
been investigated [2(| [28| less extensively as compared to 
the half- filled model. The phase diagram of the quarter- 
filled Peierls-Hubbard model is discussed in [29j |. Hence, 
we restrict our calculations to the ionic potentials. 

We access the stability regions of the ionic Hubbard 
model by introducing several potentials connected with 
different ordering processes of charge and spin degrees of 
freedom into the Hamiltonian. In particular, we inves- 
tigate which potential leads to the largest energy gain, 
and analyze against which potential or order the system 
is most unstable, indicating which order will be estab- 
lished in nature. In our calculations, we concentrate on 
the ionic Hubbard model at quarter filling, for both weak 
and strong interaction. We determine the ground state 
and the ground state energy as a function of interaction, 



doping and potential strength, using the density matrix 
renormalization group method [30(. Thus, our calcula- 
tions give detailed insights into the properties of hole- 
doped antiferromagnets. 

III. MODELS 

In the following section we first discuss the properties 
of non-interacting electrons on a chain with different pe- 
riodic potentials and the phase diagram of the homoge- 
neous Hubbard chain. In the last subsection we discuss 
some basic features of the periodic Hubbard model. 



A. Peierls model 

The Peierls model is the prototypical model to study 
the coupling of non-interacting electrons to the lattice. 
The Hamiltonian is given by 



JV 



H 



Peierls 



(1) 



where we consider a modulated hopping term with = 
t[l + u cos(Qa ■ i)]. N denotes the number of lattice sites 
and a the lattice constant. In addition, N e is the number 
of electrons on the chain and Q the wave vector of the 
periodic distortion. In case of a commensurate distortion 
we have Q = 2mfcj?, m = 1,2,.... The dimerized chain, 
Q = n/a, at half filling, n = N e /N = 1, was investigated 
by Su, Schrieffer, and Heeger (SSH) [3- Due to the 
periodic distortion a band gap, A = iut, opens at k& = 
Q/2. Since at half filling we have k& — kp, the system 
is insulating. As is characteristic for a band insulator, 
the gaps for charge and spin excitations are identical. In 
the following, we concentrate on the quarter-filled band 
and 4a-periodic distortions, because in case of the 2a- 
periodic distortion the quarter-filled system is metallic, 
since &a 7^ kp. 

To begin, we determine the energy bands and the 
energy gain due to a periodic modulation of the hop- 
ping parameter with Q = ir/2a. For quarter filling, 
kp = 7r/4a, this system shows the Peierls transition, 
since Q = ir/2a = 2k f . The Fourier transformation of 
the hopping term yields, in analogy with the SSH-model, 
the expression 



H 



E 



i- ' i'-r).» = [°k C k+Q C k+2Q C k+3Q J th a p 
k 



where the matrix h a g is given by 



C k+Q 
C k+2Q 
\ Ck+3Q 



(2) 
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/ -2 cos ka -u{c lka + ie~ ika ) u(e ika - ie~ lka ) \ 



<-a{3 



') -2sinfca -u{c~ ika + ie lka ) 



\ u(e~ lka +ie lka ) 



u(ic lka — e ) 



2 cos ka 



u(e lka + ie~ lka ) 



(3) 



u(e — ie ) 2sinfca j 



The eigenvalues of this matrix are 



e{k) 




2u 2 ± J cos 2 2ka + 4u 2 (l + sin 2 2k) 



(4) 

compare [31|. The four energy bands are separated by 
the energy gaps 



A 1 =Aut, A 2 = 2u 2 t, and A 3 = Ai, 



(5) 



where Ai = e-iik — hp) — E\{k = kp) and A 2 = £3(&: = 
2kp) — S2(k = 2kp). The energy gain of the fermions, for 
the quarter-filled case, obtained by summing up the con- 
tributions of all occupied states, is found to be cx u 2 t lnw, 
as for half filling. 



& 




B. Ionic model 

Another type of periodic distortion is given by local 
potentials, resulting in the so-called ionic model 

JV N N 

Clonic = ~tJ2 (4,v C i+l,<T+ h - cj+^WiTH-Y^hiSi . 
i.(7 i i 

(6) 

Again, for the half-filled band the same calculation as 
before can be performed for a periodic potential, with 
A c = A s = 2W . In case of the ionic model for a quarter- 
filled system we compare several potentials, which corre- 
spond to different charge and spin patterns: 

a) First, we study a simple potential with 2fcj?(4a)- 
period, given by W% = W cos(ni/2). Here, 
the charge order for zero hopping is given by 
1/0.5/0/0.5 electrons per site, periodically contin- 
ued. Spin order is not established. 

b) A modification of a) given by W{ = W cos(iri/2) + 

cos(7r(i-|-l)/2), which leads to two occupied and 
two unoccupied sites ("cluster") (l3| . i.e. a 1/1/0/0 
charge order. 

c) A charge order with a 1/0/1/0 pattern, i. e. a 
4fc^(2a)-periodic order, is forced by a potential 
with Wi = Wcos(iri). With interaction we like- 
wise expect magnetic ordering. 

d) In addition, a local magnetic field with hi = 
W 7 cos(7ri/2)(n ii | — m^i) enforces a f/0/J./O pat- 
tern. As in c), we expect a 2fcj?(4a)-period for the 
spins and correspondingly a Akp (2a)-period for the 
charges. 



FIG. 2: (Color online) One-particle energy e(k) versus mo- 
mentum k, where W = G.lt and the horizontal lines indicate 
the Fermi level of the half-filled and quarter-filled band. The 
straight line corresponds to the Q = 7r/a-periodic potential, 
the dashed line to the Q — 7r/2a-periodic potential. The re- 
duction of the Brillouin zone of the clean model is obvious. 



Diagonalizing the resulting Hamilton matrix, as de- 
scribed above, yields the dispersions 



a), d) e(k) = ±tJ2 + w 2 /2 ± ^cos 2 ^) + 2w 2 + w 4 /A, 



b) e(k) = ±tJ2 + w 2 /2± ^4cos 2 (2fc) + Aw 2 , 



c) e(k) — ±W w 2 + 4 cos 2 k, 

with w — W/t. Fig. [2]shows the dispersion for potentials 
a) and c). In comparison to the Pcierls model not only a 
gap opens at k = Q/2 but the lower and upper bands are 
shifted down and up, respectively. In case of a small per- 
turbation, the dispersion is similar for the Peierls and the 
ionic model, as well as for potentials a) and b). Potential 

a) and d) are equivalent in the non-interacting case. 
The energy gain in the quarter filled band is found to 

be about —W 165 in cases a), b), and d). For small W, 
the energy gain here is nearly the same in case a) and 

b) , whereas it becomes weaker in case b) than in case a) 
at W w 0.5i. In case c), the energy gain is quadratic, 
E(W) - E(0) = §IF 2 . It mainly traces back to the band 
shift. 

Accordingly, the gaps are given by 

a), d) Ai —W,A 2 — V2W, A 3 = W; 
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b) Ai ~ W for small W, A 1 -> Af for strong W, 
A 2 = V2-W/2, A 3 =A i; 

c) Ax = 0, A 2 = W, A 3 = 0. 

C. Hubbard model 

The Hubbard model is known to capture the interplay 
between kinetic energy (delocalization) and interaction 
(localization) in electronic systems. The Hamiltonian is 
given by 

N N 

i?Hubb = ~^2U (cf a c i+l a + h. c.) + U ^ni, T n a . 

i,a i 

(7) 

The Hubbard model in one dimension is exactly solv- 
able by means of the Bethe ansatz !l5j. Note that in 
one dimension another useful formulation of the Hub- 
bard model is available on the basis of the bosoniza- 
tion technique. The low lying excitations of the non- 
interacting as well as the interacting fermions system are 
sound waves, i. e. the Fermi system can be described as 
a non-interacting Bose system, called a Luttinger liquid, 
showing spin-charge separation. In the clean case, the 
Hubbard model has three phases. For U < 0, the spin 
excitation spectrum has a gap and the low-lying charge 
excitations can be described by those of a Luttinger liq- 
uid. For U > and away from half filling, spin and 
charge excitations are those of a Luttinger liquid. The 
last phase occurs for U > and half filling, where the 
charge excitations have a gap and the spin excitations 
are of Luttinger type. A relevant 4fc^-Umklapp scatter- 
ing term, only present for half filling, is responsible for the 
Mott gap in the latter phase. The bosonization technique 
is adequate for metallic systems or in the weak coupling 
regime. It is useful to determine the phase boundary 
between metals and insulators but it is not suitable for 
distinguishing different insulating phases for intermedi- 
ate or strong perturbations. 

D. Periodic Hubbard model 

A commensurate periodic distortion - i. e. commensu- 
rate to the band filling - introduces an additional non- 
linear term in the bosonized Hamiltonian, which couples 
spin and charge degrees of freedom and destroys the in- 
tegrability of the clean Hubbard model. In the half-filled 
case we therefore find a transition between two possi- 
ble insulating phases, i. e. the band (Peierls) insulator 
(W, ut ^> U) with charge and spin gap of the same order 
and the Mott insulator (W, ut <C U) with a charge gap 
proportional to the interaction, but vanishing spin gap, 
where periodic distortion and interaction oppose each 
other. However, the transition is fundamentally different 
for the Peierls and the ionic Hubbard model [l6|, [l7[ • In 
the Peierls model the transition between the two phases 



is smooth and continuous [18|, whereas a third phase - 
a spontaneous dimerized insulating phase - is found in 
the ionic Hubbard model. Within this phase, charge and 
spin gap are non-zero, and the dimerization operator has 
a non-zero expectation value. In fact, the extension of 
this phase is discussed controversially. 

The situation in the quarter filled band is completely 
different. With interaction, no insulating phase is estab- 
lished, except for very strong interaction, where higher 
order terms of the Umklapp scattering can cause a Mott 
insulator on its own (26j. In the quarter-filled case with 
2fcp-periodic potential we expect from bosonization a 
transition from the Luttinger liquid to the band insulator 
in the non-interacting as well as in the interacting sys- 
tem. The interaction modifies the exponents of the gaps 
and the energy gain. Since they become smaller with 
increasing interaction, interaction stabilizes the band in- 
sulator. Here, the interplay between periodic distortion 
and interaction is weak, since the quarter filled model 
with interaction is still a metal. A similar situation is 
found for half filled spinless fermions for weak interaction, 
where the Umklapp scattering induces the insulator at fi- 
nite interaction strength [27j • In case of a 4fc^-periodic 
potential Umklapp scattering becomes relevant [28| due 
to the doubling of the unit cell or, equivalently, due to 
the reduction of the Brillouin zone, see the straight line 
in Fig. 1. In this case, the system undergoes a transi- 
tion from the Luttinger liquid to a Mott insulating phase 
when the potential and the interaction are present. Po- 
tential and interaction hence cooperate rather than op- 
pose each other as in the half filled ionic model. 

An analysis of the ionic Hubbard model with another 
2fcj?-periodic potential is given in [19(. The authors find 
the band insulator for weak interaction but strong poten- 
tial. Our 4/cF-periodic pattern, hence the Mott insulator, 
is obtained in their model in the limit of strong inter- 
action. Both phases are separated by a bond ordered 
phase. Thus, the model potential used in [l9[ reveals 
the features of the half filled ionic Hubbard model. The 
analysis, however, concentrates on large \U — W\. In our 
case, the 4/cf-periodic pattern is not a limit of the 2kp- 
periodic patterns. 

In the following, we concentrate on the ionic Hubbard 
model at weak to intermediate potential strength and 
intermediate to strong interaction. In addition, we study 
the region of stability for each potential. 



IV. NUMERICAL RESULTS FOR THE IONIC 
HUBBARD MODEL AT QUARTER FILLING 

In the following we discuss the numerical data for the 
ground state energy, the charge and spin gap, and the 
spin-spin correlation function. The numerical results are 
obtained by the density matrix renormalization group 
method (DMRG), as implemented by Brune 32]. The 
DMRG is a quasi-exact numerical method to determine 
the ground state properties, i. e. the ground state and 
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the ground state energy, of long one-dimensional (non- 
integrable) systems with reasonable accuracy [3(j. Re- 
garding different boundary conditions, it is useful to take 
into account the equivalence of fermion and spin models 
and to implement the spinless-fermion model in terms of 
an equivalent spin chain, and the Hubbard model as two 
coupled spin chains without perpendicular XY-coupling. 
Using the DMRG, it is possible to extend the tractable 
system lengths for the Hubbard chain to N = 60 . . . 100 
sites. In our simulations we perform five lattice sweeps 
and keep 300 to 500 states per block. Correlation func- 
tions can be obtained within an error of 10 -6 in the Hub- 
bard model when using open boundary conditions. A 
memory of about 700 MB is required. 

A. Energy 



b --X- 
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First, we determine the ground state energy as a func- 
tion of W and U. The energy gain depends almost 
quadratically on the potential strength for all potential 
types, and is dominated by the band shift. The numerical 
data are depicted in Fig. [3l In case of a coupling to the 
charge density, Willi, the 2/cF-periodic potential is stabi- 
lized for small U (case a)), and the 4/cF-periodic potential 
(case c )) f or large U, as found for the Hubbard-Holstein 
model [331 ] . This energy gain is mainly due to the prefac- 
tors. A more detailed analysis of the algebraic behavior 
for W ->■ shows that E a cx FT 168 and E c cx W 2 if 
U = 2t, but E a cx FT 176 and E c cx W 1 - 76 if U = 10*. In 
case a) the exponent increases with interaction, in case 
c) it decreases. For all interactions, however, term d), 
hiSf , is dominant. Of course, in a real system, a poten- 
tial coupled to the density is much bigger than a magnetic 
field. It leads - like potential c) - to 4/cF-oscillations in 
the charge density and 2/cF-oscillations in the magneti- 
zation. In the non-interacting system the data points for 
potentials a) and d) lie nearly on top of each other (by 
use of open boundary conditions the exact equivalence of 
a) and d) is waived), but in the interacting system the 
energy gain due to potential d) grows faster with inter- 
action than in all other cases. To be more specific, we 
find E d cx W 1 - 62 for U = 2t and E d cx W 1A4 for U = 6i. 
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FIG. 3: (Color online) Energy gain per site versus square of 
the potential strength (both in units of t). We compare po- 
tentials a), b), c), and d). The number of sites varies between 
N = 28 and N = 60, where [E(0) - E{W)]/N is almost in- 
dependent of N. The interaction is U = It in the upper plot 
and U — 6t in the lower plot. 



B. Charge and spin gap 

We note that the excitation gaps 

A c = ~[£(iV e + 1) + E(N e - 1) - 2E(N e )} 

and 

A s = ^[E(S Z = 1) + E(S Z = -1) - 2E(S Z = 0)] 

are calculated for finite systems only. An extrapolation to 
the thermodynamic limit, A = A +f(N) with f(N) — > 
for N — > co, is not performed, for the following reasons: 



for small gap Ao the accuracy of the DMRG away from 
half filling is not accurate enough and for large Ao we 
have f(N) as for the considered system sizes. The de- 
pendence of the excitation gaps on the potential strength 
is shown in Fig. [5] In cases a), b) and d) we see that 
A s ~ A c , according to the band insulating state. An 
almost linear dependence on W, as for non-interacting 
particles, is recovered in the interacting system for the 
charge gap in cases a), b), and c). The saturation of the 
gaps in case b), as discussed in Sec. 3B, is obtained for 
the spin gap in the considered parameter regime. The 
influence of the interaction in the band insulating state 
can be summarized as follows. In case a) the charge 
gap becomes smaller with interaction, the spin gap even 
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FIG. 4: (Color online) Charge and spin gap versus potential 
strength (both in units of t). In the upper plot we show A c 
and in the lower A a , comparing potentials a), b), c), and d). 
The number of sites varies between TV = 28 and TV = 60 and 
the interaction strength is U = 6t. 



C. Correlation Function 

In order to calculate the spin-spin correlation function 
(S*Sj) within DMRG we again use open boundary con- 
ditions and the Parzen filter function [3iJ to reduce the 
Friedel oscillations. In the Mott insulating regime, the 
Friedel oscillations in the spin sector are long ranged, 
but decay fast in spin gap systems [35j |. The main ques- 
tion concerns the correspondence of the spin-spin corre- 
lation functions of the 2a-periodic potential (case c)) and 
the Heisenberg model. For comparison, we calculate the 
spin-spin correlation function of the half-filled Hubbard 
model, which can be mapped onto a Heisenberg chain in 
case of strong interaction. In the Heisenberg model we 
have 



(S i S ) — 



1 



o-2k F cos(2kpa ■ i) 



47r(a • i) 2 



In gapped systems, the spin-spin correlation functions 
reflect the behavior of the energy gaps. Thus, the spin- 
spin correlation function decreases exponentially in cases 
a) and b), which show a spin gap without spin order. 
In case d) the spins are fixed by the potential hi. The 
exponentially decaying part is hard to extract from the 
mean magnetization in the numerical data due to the 
incomplete suppression of the boundary oscillations. 

In the Mott insulating regime, however, the linear de- 
crease of the spin-spin correlation function with distance 
is obtained already for small potential W. For bet- 
ter comparison with the data of the half-filled Hubbard 
chain, we show data for W — t in Fig. [51 Obviously, 
the quarter-filled Hubbard chain with 4kp-periodic dis- 
tortion can be mapped onto the Heisenberg model, too. 
We note that in case d) already a weak potential - con- 
nected with a small spin gap - leads to significant struc- 
tures in (S*Sj), whereas in case c) the potential has to 
be much stronger to yield effects of similar size. The 
increase of S(q) oc q appears both in the half-filled Hub- 
bard chain and in the quarter-filled Hubbard chain with 
4fcp-periodic potential. S(q) has a sharp maximum at 
q = ir/2 in case d). 



more. In case b) the charge gap increases with inter- 
action for small U, but saturates for strong interaction. 
The spin gap decreases rapidly with interaction. In case 
d) the charge gap increases but does not saturate. The 
spin gap shows a maximum for intermediate interaction 
strength. For U = lOt it is even smaller than in the 
non-interacting system. On the contrary, we find for po- 
tential c) a linear increase of A c with U and W as well as 
A S (W, U, TV) - l/U, 1/W,1/N -> 0, indicating the Mott 
insulating state. Due to the finite system size, we are not 
able to obtain exponents for small IF. As a consequence, 
a comparison with bosonization or the energy data, see 
Section A, is not possible. To conclude, for strong inter- 
action and small potential we find a large charge gap but 
a small spin gap. 



V. SUMMARY 

In summary, we have studied the effects of periodic 
potentials on a chain. In particular, we have considered 
the one-dimensional Hubbard model at quarter filling. 
We have compared periodic potentials yielding different 
types of behavior: a) a 2fcj?-periodic potential leading 
to a band insulator, b) a 2fcj?-periodic potential leading 
to a band insulator with cluster-like arrangements of the 
charges, c) a Akp -periodic potential leading to a Mott 
insulator with antiferromagnetic alignment of Heisenberg 
type spins on next-nearest neighbor sites, and d) a 2kp- 
periodic magnetic field leading to a band insulator with 
antiferromagnetic alignment of Ising type spins. Cases c) 
and d) reveal the same charge and spin distribution but 
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FIG. 5: (Color online) Spin-spin correlation function {S* S34} 
versus site i. We compare the half- filled Hubbard chain 
(straight line; which is equivalent to a Heisenberg chain), 
with the quarter-filled Hubbard chain on a 4fcF-periodic po- 
tential (long-dashed line, W = t; equivalent to a half-filled 
Heisenberg chain), and the quarter- filled Hubbard chain in 
a 2fcF-periodic magnetic field (short-dashed line, W = O.lt; 
equivalent to a half- filled Ising chain). In all data sets the 
interaction is U — 5t. 

different spin excitations. In the Heisenberg case the spin 
excitations are gapless, whereas in the Ising case they are 
gapped. 

In non-interacting systems, we find E a — Ed > Eb > 
E c , thus the band insulator. The charge distribution is 
quite homogeneous, and the cluster like arrangement is 
not minimal in energy. Turning on the interaction, poten- 
tial d) results in the largest energy gain, while the order of 
the remaining potentials depends on both the potential 
strength and the interaction. The interaction strongly 
supports the spin order, where double occupancy is sup- 
pressed. 

In the following discussion we rely on the potentials a) 
to c), where the spin-order, if present, is a result of the 
charge order. For weak interaction and weak potential, 
we recover the behavior of the non-interacting system, 
with E a > Eb > E c . Turning to intermediate values 
(U 1=3 3t, W 1=3 t ), the cluster-like arrangement of poten- 



tial b) gains more energy than the homogeneous distri- 
bution of case a,), Eb > E a > E c . For strong W — ► t, also 
a 4kp pattern is favored against the homogeneous case, 
Eb > E c > E a . Thus, for small or intermediate interac- 
tion and potential, the gain of energy due to the hopping 
is stronger than the effects of the repulsive interaction 
or the underlying potential. With increasing potential, 
however, the effects of the potential dominates. 

On the other hand, regarding a weak potential but in- 
creasing interaction, we find that a pattern related to 
potential a) is established only at small interaction. The 
cluster-like arrangement is formed at intermediate inter- 
action (U « 3 — 5i), but at strong interaction the 4fcp- 
periodic pattern dominates. Concerning the spin chain 
compound (Ca,Sr)i4Cu2404i, we have applied the po- 
tential only in order to clarify the leading instability of 
the system. For comparison with a real material which 
shows no lattice effects, the limit W — > is relevant. 
In addition, for copper oxides usually a correlation pa- 
rameter of U « 8eV [36| is assumed. Together with the 
band width of about 1 eV, see Fig. 1, we have U > I5t. 
In this parameter region, we find - neglecting the mag- 
netic field - the 4fcp-periodic charge pattern with anti- 
ferromagnetism to be the most probable ground state. 
The energetic order found for weak potentials extends to 
strong potentials. For strong interaction, we always have 
E c > Eb > E a . The tendency of the repulsive interaction 
to separate the charges dominates the phase diagram in 
this region. The influence of the potential is weaker, lead- 
ing to a preference of the cluster over the homogeneous 
distribution. 

Finally, we remark that adding electrons stabilizes the 
cluster formation. This observation likewise agrees with 
experimental results, obtained for the (Ca,La)i4Cu2404i 
series, where the cluster size grows on additional charge. 
Therefore we conclude, that the antiferromagnetism is 
due to the band filling, thus a commensurability effect. 
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